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^ , ABSTRACT 

C \ We present a detailed analysis of the complex soft X-ray spectrum of NGC 4151 

~^ ' measured by the RGS instruments aboard XMM-Newton. The XMM-Newton RGS 

spectra demonstrate that the soft X-ray emission is extremely rich in X-ray emission 
I lines and radiative recombination continua (RRC), with no clear evidence for any 

underlying continuum emission. Line emission, and the associated RRC, are clearly 
detected from hydrogen-like and helium-like ionization states of neon, oxygen, nitrogen 
and carbon. The measured lines are blueshifted with a velocity of between ^100-1000 
' km s~^, with respect to the systemic velocity of NGC 4151, approximately consistent 

ly-^ I with the outflow velocities of the absorption lines observed in the UV spectrum of NGC 

l/^ , 4151 (Kriss et al. 1995), suggestive of an origin for the UV and soft X-ray emission in 

the same material. Plasma diagnostics from the observed helium-like triplets, imply a 
I range of electron temperatures of '^1-5x10* K and electron densities of between 10*- 

'nJ" . 10^° cm~^. The soft X-ray spectrum of NGC 4151 is extremely similar to that of NGC 

' 1068, both in terms of the atomic species present and in terms of the relative strengths 

fH , of the observed emission lines and RRC (Kinkhabwala et al. 2002), suggesting that the 

Q^' soft X-ray excesses observed in many Seyfert 2 galaxies may be composed of similar 

emission features. Modelling the RGS spectra in terms of emission from photoionized 
and photoexcited gas in an ionization cone reproduces all of the hydrogen-like and 
helium-like emission features observed in the soft X-ray spectrum of NGC 4151 in 
^ I detail and confirms the correspondence between the soft X-ray emission in NGC 4151 

and NGC 1068. NGC 4151 shows somewhat lower individual ionic column densities 
than those measured for NGC 1068 indicating that the material in the ionization cones 
of NGC 4151 may be somewhat more dense than the material in the ionization cones 
of NGC 1068. 
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1 INTRODUCTION 

Prior to the advent of XMM-Newton and Chandra, the ori- 
gin of the soft X-ray emission from heavily obscured AGN 
remained very uncertain. Recently, however, observations 
with XMM-Newton and Chandra have revolutionized our 
understanding of heavily obscured AGN (for a complete de- 
scription of XMM-Newton, Chandra and their relative ca- 
pabilities see Jansen et al. 2001, Weisskopf et al. 2002 and 
references therein). The combination of the high spatial res- 
olution afforded by Chandra, the large collecting area of 
XMM-Newton and the high spectral resolution X-ray spec- 
troscopy capabilities (in the form of X-ray grating and CCD 
instruments) of both missions across the 0.1-12 keV band 
is beginning reveal many new details of the X-ray emis- 



sion from heavily obscured AGN, particularly with regard 
to the soft X-ray emission from these objects. Chandra ob- 
servations of several heavily obscured AGN have revealed 
considerable extended soft X-ray emission in each source, 
that is tightly correlated with the O[III](A5007) emission 
observed with the Hubble Space Telescope (HST; NGC 1068 

- Young et al. 2001; Mrk 3 - Sako et al. 2000; NGC 4151 - 
Ogle et al. 2000; Circinus - Sambruna et al. 2001). In sev- 
eral of these sources {e.g. NGC 1068 - Young et al. 2001), 
the extended X-ray emission and the O[III](A5007) emission 
clearly traces the spiral arms and other large-scale structures 
out to a considerable distance from the active nucleus of the 
galaxy. Similarly recent spectral observations of these AGN 
with the X-ray grating instruments on Chandra (NGC 1068 

- Brinkman et al. 2002, Ogle et al. 2002; Mrk 3 - Sako et 
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al. 2000; NGC 4151 - Ogle et al. 2000; Circinus - Sambruna 
et al. 2001) and XMM-Newton (Kinkhabwala et al. 2002) 
have revealed soft X-ray spectra dominated by prominent 
emission lines from a range of light elements {e.g. C, N, O, 
Ne etc). These observations constitute a growing body of 
evidence that suggests that the soft X-ray emission from 
heavily obscured Seyferts is largely composed of soft X-ray 
emission lines that originate in gas that is photoionized (and 
photoexcited) by the central AGN. 

The Seyfert 1.5 galaxy NGC 4151 was identified as an 
X-ray source over thirty years ago (Gursky ct al. 1971). One 
of the brightest AGN accessible in the X-ray band, NGC 
4151 has been extensively studied by missions such as EX- 
OSAT, Ginga, ROSAT, ASCA, CGRO and more recently 
RXTE, BeppoSAX, Chandra and XMM-Newton. This ob- 
servational focus has revealed that the spectrum of NGC 
4151 from 0.1-100 keV is comprised of a complex mixture 
of emission and absorption components, probably originat- 
ing in a variety of locations from the innermost parts of 
the putative accretion disk in NGC 4151, out to the ex- 
tended narrow-line region of the galaxy. Despite the rel- 
atively poor quality of the spectra, previous observations 
(particularly observations with ASCA) did tentatively re- 
veal several unresolved soft X-ray features, suggesting that 
the soft X-ray emission was somewhat more complex than 
had been previously suggested. The most prominent feature 
found with ASCA was at an energy of ~0.88 keV and was 
tentatively associated with O VIII Lya RRC (Weaver et 
al. 1994; Griffiths et al. 1998). The detailed nature of the 
soft X-ray emission in NGC 4151 has only recently been 
revealed through a Chandra HETG observation (Ogle et 
al. 2000, Yang et al. 2001). This observation demonstrated 
that the soft X-ray spectrum was composed almost entirely 
of emission linos and RRC, with little or no detected con- 
tinuum emission of any form (and hence no absorption fea- 
tures). Utilizing the excellent spatial resolution of Chandra, 
Ogle et al. (2000) resolved ~70% of the soft X-ray emis- 
sion as extended emission, originating up to 11.5" from the 
AGN, that demonstrated an excellent correspondence with 
the bi-conical O[III](A5007) emission previously imaged with 
the HST (Evans et al. 1993). The Chandra observation was 
even capable of resolving the spatial profiles of several of 
the brightest emission lines (including the iron Ka emission 
line) showing that the individual line emission was genuinely 
extended along the direction of the optical ionization cones. 
The Chandra observation dramatically confirmed that the 
hard X-ray spectrum of NGC 4151 is strongly cut-off below 
2 keV and that the soft excess is not predominantly pro- 
duced by electron scattering of the hard X-ray continuum. 
The Chandra analysis finds evidence for both photoionized 
and coUisionally ionized (thermal) gas in the ENLR, however 
many of the details of the ionized medium remained elusive, 
due to the large number (33 in total) and low significance of 
the features detected. 

This paper is organized as follows. Section 2 details the 
XMM-Newton observations of NGC 4151. Section 3 presents 
the XMM-Newton RGS spectra and discusses the spectral 
modelling in terms of emission from an ionization cone. Sec- 
tion 4 compares the RGS spectra of NGC 4151 with the soft 
X-ray spectra of NGC 1068 and Mrk 3. Section 5 presents a 
preliminary look at the correspondence between the XMM- 



Newton RGS and the XMM-Newton EPIC soft X-ray spec- 
tra of NGC 4151. 



2 THE XMM-Newton OBSERVATIONS 

NGC 4151 was observed with XMM-Newton five times dur- 
ing December 2r*-23'"'*, 2000. The RGS instruments were 
active for all five observations and the data were processed 
with the standard RGS pipeline processing chains incorpo- 
rated in the XMM-Newton SAS v5.4.1. Dispersed source 
spectra and background spectra were extracted with the au- 
tomated RGS extraction tasks in the SAS. The SAS filters 
the RGS events in dispersion-channel vs CCD-pulse-height 
space to separate the data from the different spectral orders. 
Investigation of the individual spectra indicate that the 2"'' 
order RGS spectra contain a negligible amount of informa- 
tion (a count rate of <10% of that in the 1"* order spectra, 
despite having ~60% of the effective area below lOA) and as 
a result only the 1^* order spectra are presented here. The re- 
sulting spectra have a total effective exposure time of ~128 
ks and ~125 ks for the RGS 1 and 2 instruments respec- 
tivly. A previous investigation of the XMM-Newton EPIC 
spectra confirms that the soft X-ray spectrum of NGC 4151 
remained relatively constant for the duration of the obser- 
vation (Schurch et al. 2003). The spectra from the five ob- 
servations were co-added to produce a single high signal-to- 
noise source spectrum and a single background spectrum for 
each RGS instrument. The full spectra and individual linos 
were analyzed with both maximum likelihood and min- 
imization techniques. For the maximum likelihood analysis 
the spectra were left unbinned however, for the mini- 
mization techniques (which, in general, are not particularly 
well-suited to fitting the sharp features observed in this and 
similar high resolution soft X-ray spectra) the spectra were 
binned to a minimum of 20 counts per spectral channel. 
Investigation of the spatial extent of the X-ray emission ob- 
served by XMM-Newton reveals that the soft X-ray emis- 
sion is not significantly extended on scales greater than the 
mirror PSF, consistent with the details of the previously re- 
ported Chandra observation of NGC 4151 (Ogle et al. 2000). 



3 MODELLING THE RGS SPECTRA 

Figure 1 shows the XMM-Newton RGS spectra from the 
combined observations of NGC 4151. The spectrum is com- 
pletely dominated by line emission and radiative recombina- 
tion continua (RRC). Line emission is clearly detected from 
hydrogen-like and helium-like ionization states of neon, oxy- 
gen, nitrogen and carbon and in many cases the line emis- 
sion associated with these ions is not limited to the Lya 
(2p-ls) transition but includes higher order resonance tran- 
sitions (np-ls) labeled /3-e for the different hydrogen- like and 
helium-like species respectively. Much of the unidentified 
line emission below ~16A can be associated with a mul- 
titude of unresolved iron L-shell emission lines, along with 
some contribution from hydrogcii-like and helium-like silicon 
and magnesium. The spectra show no significant continuum 
emission and no fluorescent line emission (there may be a 
marginal detection of silicon Ka fluorescence at ~7.lA, how- 
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Figure 1. RGS 1 (red) and RGS 2 (blue) spectra of NGC 4151. The spectra here are unbinned. Spectral discontinuities are caused by chip 
gaps in the CCD arrays and the effects of the damaged RGS CCDs. Unambiguously identified hydrogen-like and helium-like emission lines 
are labeled along with the clearly identified RRC. 



ever without a good knowledge of the iron L-shell emission 
in this region it is not possible to measure the line cleanly). 

3.1 The emission lines 

The spectrum is remarkably similar to the soft X-ray spec- 
trum of the archetypal bright Seyfert 2 galaxy, NGC 1068, 
reported by Kinkhabwala et al. (2002). Both sources show 
strong emission line triplets from helium-like nitrogen and 
oxygen, with a similar pattern of resonance, intercombina- 
tion and forbidden line strengths. Following the analysis 
of the XMM-Newton RGS observation of NGC 1068 pre- 
sented in Kinkhabwala et al. (2002), line fluxes, centroid 
energies and line widths were determined for the bright, un- 
ambiguously identified, emission lines observed in the spec- 
tra. The spectral fits quoted here for individual strong line 
features were performed with minimization techniques 
(and correspondingly the data were binned to a minimum 
of 20 counts per spectral channel). An examination of the 
fit properties using maximum likelihood fitting techniques 
for individual lines reveals that the best-fit parameters are 
not significantly different using these statistical techniques 
in comparison with the fits, as expected given the bright 
nature of the features examined here. The intrinsic profiles 
of the emission lines are represented by a line model with 
a Gaussian profile (we caution, however, that the original 



line profile may in fact not be Gaussian, although it does 
appear to be a good approximation for the lines presented 
here) which allows the determination of accurate line cen- 
troid energies, fiuxes and, in some cases, line widths. The 
results of fitting this profile to the brightest emission lines 
are given in Table 1. 

The measured fiuxes take account of the Galactic col- 
umn density (2.1x10^" cm"^; Warwick et al. 1995) using 
the tbabs model in XSPEC, and the line shifts are mea- 
sured with respect to the source redshift of z=0. 00332. The 
systematic uncertainty on measuring the centroid line ener- 
gies of the lines is estimated to be ~8xlO~'^A (which corre- 
sponds to about 150 km s~^; XMM Users Handbook v2.1, 
2002). The quoted errors are 90% confidence levels as de- 
fined by a Ax^=2.71 criterion {i.e. assuming one interesting 
parameter) and include the systematic uncertainty added 
in quadrature. Most of the measured line centroid energies 
are consistent with a mild blueshift, of between ~100-1000 
km s~^, with respect to the systemic velocity of the lines. 
This implies a line-of-sight outflow velocity for the ionized 
material of ^1000 km s~^, consistent with the line shifts of 
the absorption lines observed in the UV spectrum of NGC 
4151 (~600 km s"^ Kriss et al. 1995). Many of the ob- 
served line widths are largely unconstrained, with the spec- 
tral fitting only yielding upper limits (these are not quoted 
in Table 1); however the seven line widths that are con- 



© 0000 RAS, MNRAS 000, 000-000 



4 N.J. Schurch, R.S. Warwick, R.E. Griffiths & S.M. Kahn 

Table 1. Line properties of individual emission lines observed in the RGS spectra of NGC 4151 



Line ID 


^Obs 

A 


\a 


Line Shift 
km s^i 


"'Obs 

xio-2 A 


Flux 

Xl0~4 Photons s~i cm~2 


Ne X a 
Ne IX r 
Ne IX i 

N(- IX r 


12 13+0-01 

, o 40+0.02 

10 co+0.01 
lO.0O_Q 02 


12.134 
13.447 
13.550 

13.697 


^^"-250 
.430+490 

-470+24° 

^'"-340 

-310+22" 
-'1^' .L8IJ 


2.6li;^ 

2.411;? 
2 0+1-° 

Z.U_j 2 


42+0-05 
47+0.05 

"■^'-0.05 

46+0-1^ 

1 03+0-1 


O VIII p 
O VIII a 
O VII r 
O VII i 
O VII / 


16.00100^ 
18.96«:°1 
21 59+0-01 

21.79«;«1 
22.081H1 


16.006 
18.969 
21.603 
21.798 
22.101 


-150+230 
210+130 
230+"° 
-230l}f„ 
-300111° 


1+0.6 
0.68l°;°«9 


07+0.04 

1 60+0-10 

1 19+0.09 

QQ + 0.17 

•3.09_o.l9 


N VII 13 
N VII a 
N VI r 
N VI i 
N VI/ 


20.891°;°! 
24.77t0:°l 

28.751°;°} 
29 09+0-°^ 

29.511°°} 


20.910 
24.781 
28.792 
29.075 
29.531 


-230+1°° 

-190111° 
-350111° 

qn+250 
•^^-160 

-290+1°° 


0.71°;^ 


22+°-°* 

r, O1+0.09 

0-»i-o.os 
45+0-07 

S6+0-°° 

1 i';+o-i2 


C VI /3 
C VI a 


28.4310;°} 
33.691°;°} 


28.466 
33.736 


oQO+130 


0.9l°;| 


n 07+0.05 
"■"^'-0.05 
, 70+0.16 



" Vainshtein & Safronova, 1978 



strained by the spectral fitting appear to show a bimodal 
distribution -with the shorter -wavelength lines having an av- 
erage line width of o-~0.025 A(~1000 km s^i FWHM) and 
the longer wavelength lines having an average line width of 
~0.005 A(~200 km s"^ FWHM). The observed Une widths 
imply that the intrinsic velocity dispersion of the outflow- 
ing material is <^1000 km s~^ which is somewhat lower than 
the lino widths of the absorption linos observed in the UV 
(1000-2500 km s-^ Kriss et al. 1995). The bimodal distri- 
bution of the line widths as a function of the wavelength of 
the line is suggestive of a scenario in which individual clouds 
of emitting material close to the source of the ionizing con- 
tinuum are moving more rapidly than clouds further from 
the continuum source. The large majority of the soft X-ray 
flux of NGC 4151 as measured by Chandra is not extended 
on scales larger than ~2", demonstrating that the reported 
extended soft X-ray flux will not contribute considerably to 
the line widths measured in the RGS spectrum. 

3.2 Temperature diagnostics -with Radiative 
Recombination Continua 

The radiative recombination continua (RRC) identifled in 
the RGS spectra are the result of a free electron recombin- 
ing directly into a bound state of an ion. Recombination to 
the ground state represents the dominant transition and is 
the only transition considered during the spectral modelling. 
The free electrons can have any energy above the recombi- 
nation threshold and hence the feature produced is neither 
symmetrical or narrow. The energy of the electrons prior to 
the recombination process is dependent upon the electron 
temperature of the emitting plasma and hence the width 



of the RRC features provides a direct measurement of the 
recombining electron temperature, Tc (Liedahl & Paerels, 
1996; Liedahl 1999). RRC originating in hot plasmas are 
broad and often subtle features however, RRC from cool 
plasmas are narrow, prominent, features. The velocity dis- 
persion of the emitting material will affect the RRC in ex- 
actly the same manner as it affects the emission lines dis- 
cussed previously, however this effect is negligible here due 
to the low velocities of the ions involved compared to that 
of the electrons. RRC from hydrogen-like and helium-like 
ionization states of No, O, N and C are clearly identifiable 
in Figure 1. Electron temperatures were determined for the 
bright, unambiguously identified, RRC (Table 2) using the 
RRC profiles generated by the 'ionization cone' model dis- 
cussed in Section 3.4 fit, using minimization techniques, 
to the individual RRC. Again, an examination of the fit 
properties using maximum likelihood fitting techniques for 
individual features reveals that the best-fit parameters are 
not significantly different using these statistical techniques 
in comparison with the fits. The derived electron tem- 
peratures imply electron temperatures of ~2-4xl0'' K. The 
association of high- ionization emission lines {e.g. O VIII, 
Ne X etc; identified in the previous section) with such a low 
temperature plasma is strongly indicative of a plasma whose 
ionization state is governed by photoionization as opposed 
to coUisional ionization. 

3.3 Plasma diagnostics with helium-like triplets 

The three most prominent X-ray emission lines from helium- 
like ions correspond to transitions from the n=2 shell to the 
n=l shell and are the resonance (ls2p ^Pi^ls^ ^So), in- 
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Table 2. Electron temperatures (in eV) from RRC 



RRC 7e 



C V 


^■^-0.5 


C VI 




N VI 


O Q + 3.0 


N VII 


4 0+2-° 

^•"-1.5 


O VII 




O Mil 


2.0 



Table 3. Helium-like triplet plasma diagnostic ratios 



Triplet ID 


R 


G 


NelX 


9 Q+0.7 
^■•^-0.7 


Q 9+1.2 
■^■^-0.9 


OVII 






NVI 




3 4+1.0 



tercombination (ls2p ^P2,i— >ls^ '^So) and forbidden (ls2s 
^Si^ls^ ^So) line transitions (Pradhan ct al. 1981). Ratios 
of the measured line fluxes for the recombination, intercom- 
bination and forbidden lines from helium-like ions provide 
a diagnostic of the electron density and temperature in hot 
plasmas. Gabriel & Jordan (1969) were the first to recog- 
nise the importance of the helium-like triplet lines and used 
line strength ratios to analyze the solar corona spectra in 
detail. The ratio of the intensity of the forbidden lino to the 
intensity of the intercombination line (R, where R=f/i) is 
sensitive to the electron density of the plasma and the ratio 
of the sum of the intensities of the forbidden and intercombi- 
nation lines to the intensity of tfie resonance line [G, where 
G=(i-|-f)/r] is sensitive to the electron temperature of the 
plasma. For a detailed discussion of the manner in which 
the ratios depend on the electron density and the plasma 
temperature see Pradhan et al. (1981), Porquct & Dubau 
(2000), Porquet et al. (2002) and references therein. 

The R and G ratios for the Ne, O and N triplets ob- 
served in the RGS spectrum of NGC 4151 are given in Ta- 
ble 3. The G-values imply an electron temperature of <W^ 
K, consistent with the electron temperatures implied from 
the properties of the observed RRC and consistent with pho- 
toionization being the dominant ionization mechanism in 
the plasma (Porquet & Dubau, 2000). The R- values imply 
electron densities of ~10*-10'^ cm^'^ in nitrogen and oxygen 
and ~10^°-10^^ cm"-^ from neon (Porquet & Dubau, 2000). 
In particular the ratios calculated from the oxygen triplet 
(which is both the strongest and the leeist confused of the 
helium-like triplets) strongly imply a photoionizod plasma 
however, extracting detailed values of density and tempera- 
ture for the emitting material is complicated by the possi- 
bility that a strong radiation field {i.e. photocxcitation) can 
mimic the effect that a relatively high density has on the 
line strengths. In such a situation the electron density, and 
to a lesser extent, the electron temperature can be overes- 
timated (Porquet et al. 2002). We note that the expected 
contribution of blended dielectronic satellite lines to these 



Seyfert 1 



, {=n/4n 




Figure 2. Schematic cartoon of the 'ionization cone' model 
(not to scale) . The central nuclear components (black hole and 

accretion disk), denoted by the black spot, arc surrounded by 
a dense absorbing torus of neutral dust and gas. The optically 
thin material perpendicular to the plane of the absorbing torus 
sees the direct continuum from the X-ray source and becomes 
strongly ionized forming ionization cones. The line emission 
from these cones is observable from Seyfert 2 galaxies due 
to the suppression of the direct continuum, whereas the ion- 
ization cones in Seyfert 1 galaxies axe only seen through their 
absorption of the direct continuum (Kinkhabwala et al. 2002). 

ratios has not been taken into account in the values quoted 
here. For the N and O triplets this contribution is negligible; 
however the Ne triplet is significantly affected. For a high 
temperature, collisionally dominated, plasma the additional 
contribution can result in an error of ^1% and ~9% for 
the R and G ratios respectively, however for a photoionized 
plasma, where the temperatures are several orders of magni- 
tude lower, this effect can be considerably greater (Porquet 
et al. 2002). 

3.4 The ionization cone model 

Given the similarity between the RGS spectra of NGC 1068 
and NGC 4151, wc utilize the 'ionization cone' model devel- 
oped by Kinkhabwala et al. (2002) in their analysis of the 
RGS spectrum of NGC 1068, to fit the RGS spectrum of 
NGC 4151. This model, which was incorporated into XSPEC 
for ease of spectral fitting, reproduces the hydrogen-like and 
helium-like line emission, and the associated RRC, from an 
irradiated cone of optically thin plasma (c./. the ionization 
cones observed in 0[III] by HST) shown schematically in 
Figure 2. The model includes the effects of both photoion- 
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ization and, crucially, photoexcitation. Photoexcitation is 
the process by which the absorption of a photon with insuffi- 
cient energy to ionize the atom, results in the transition of an 
inner-shell electron to a higher energy level within the atom. 
The excited electron will rapidly undergo one or more radia- 
tive transitions, decaying back to its original energy state. 
For each atomic species, photoexcitation is expected to en- 
hance the intensity of the higher order resonance transition 
lines relative to the order resonance line. Photoexcita- 
tion was observed to play a significant role in the spectrum 
of NGC 1068 where the Ly/3, Ly7 and LyJ lines were sig- 
nificantly stronger than predicted by models based on pure 
photoionization or coUisional ionization. 

Here, as with NGC 1068, the plasma is irradiated by 
the unabsorbed intrinsic power-law continuum, the slope of 
which in the case of NGC 4151 is fixed at r=1.65, based on 
the analysis of the hard X-ray XMM-Newton EPIC spec- 
trum presented in Schurch et al. (2003). The luminosity 
of the central source, Lx, and the covering fraction of the 
cone, f, are fixed at Lx=6xl0^'' erg s^^ and f=0.1, respec- 
tively. The modelling includes both Galactic column den- 
sity (2.1x10^° cm"^; Warwick et al. 1995) using the tbabs 
model in XSPEC and a source redshift of z=0. 00332. The 
electron temperatures for each atomic species were fixed at 
the temperatures found from the analysis of the RRC (Ta- 
ble 2), where possible, and the remaining RRC temperatures 
were fixed at the temperatures identified for NGC 1068 by 
Kinkhabwala et al. (based on the similarity between the ob- 
served bright RRC in NGC 4151 and the corresponding RRC 
in NGC 1068). The remaining free parameters of the model 
are the transverse velocity distribution of the lines, cr"''^, the 
radial velocity distribution of the lines, crj"'*, and the radial 
column density of each individual atomic species, N™^. The 
modelling presented in the next section includes a significant 
contribution from iron L-shell emission lines. Iron L-shell 
emission, known to be present in the soft X-ray spectrum 
of both Mrk 3 and NGC 1068, will contribute significantly 
between ~9-16A(0.75-1.4 keV) and will results in an arti- 
ficially high measure of the ionic column densities of Neon 
and Magnesium if it is not incorporated in the modelling. 

3.5 Fitting the ionization cone model to the RGS 

spectra 

The ionization cone model was fitted to the RGS spectra 
in XSPEC for simplicity. We note that the minimization 
techniques, used earlier in the paper to fit specific bright 
spectral features, are, in general, not well-suited to fitting 
wide-band X-ray grating spectra which often incorporate 
many additional low-significance spectral features. With this 
in mind, the ionization cone model is fit to the RGS spectra 
of NGC 4151 using maximum likelihood fitting techniques. 
The complexity of this specific model makes the accurate 
determination of individual parameter errors problematic. 
Despite this caveat, the model provides a "good fit" to all 
the hydrogen-like and helium-like line emission and RRC 
observed in the RGS spectra. Here, the term "good fit" here 
refers only to a set of parameters that correctly reproduces 
the distinct features observed in the spectrum (i.e. it is not 
a statistical measure of the goodness of fit of the model). 
Because of this, statistical errors are not quoted for either 
the column densities or the velocity measurements. Figure 3 



Table 4. Best-fit parameters for the 'ionization cone' model 
fit to the RGS spectra of NGC 4151 



ion 


N™'' (cm-2) 

ton ^ ' 


kTe (eV) 


C V 


2.9x10" 


2.5 


C VI 


4.0x10^^ 


4.0 


N VI 


1.4x10^^ 


3.0 


N VII 


2.4x10" 


4.0 


O VII 


7.6x10 ' 


4.0 


O VIII 


7.9x10" 


6.0 


Ne IX 


4.6x10" 


4.0" 


Ne X 


4.7x10" 


4.0" 


Mg XI 


2.2x10" 


4.0" 


Mg XII 


4.0x10" 


4.0" 


Si XIII 


6.2x10" 


4.0" 


Si XIV 


l.lxlOi* 


4.0" 



" Unconstrained parameter. Value 
fixed at the value used in 
Kinkhabwala et al. (2002). 



shows the best-fit of the ionization cone model to all the 
hydrogen-like and helium-like emission lines and RRC in 
the RGS spectrum of NGC 4151. The best-fit radial ionic 
column densities and the associated RRC electron tempera- 
tures are given in Table 4; a"*"" has a best-fit value of ^^^400 
km s~^ and cr™'* has a best-fit value of --^330 km s~^. We 
do not quote any specific best-fit values for the ionic column 
densities of iron responsible for the iron L-shell emission (Fe 
XVI- XXIII) due to the unresolved nature of the specific iron 
L-shell lines. Typically the ionic iron column densities are of 
the order of ~few xlO^® cm~^. 

Both the velocity measurements axe in reasonable agree- 
ment with the average line widths and line shifts measured 
from the individual line centroid energies. Kinkhabwala et 
al. (2002) quote, based on the NGC 1068 data, an estimate 
for the uncertainties on the individual model parameters of 
approximately a few tens of percents for the ionic column 
densities and approximately a factor of 2 for the transverse 
and radial velocities. Here, the data are of somewhat worse 
statistical quality than the NGC 1068 spectra and so we ex- 
pect the parameters to be less well constrained. NGC 1068 
is a factor of ~3 times brighter than NGC 4151 and hence 
the errors on the fit parameters from the NGC 4151 spec- 
trum will be a factor of ~2 times greater than the errors 
estimated for the NGC 1068 fit parameters. Using a frac- 
tional ionic abundance of 0.5 for O VIII, the column density 
of O VII quoted in Table 4 yields an equivalent hydrogen 
column density of Njf^2xl0^^ cm~^. 



4 A COMPARISON WITH NGC 1068 AND 
MRK 3 

The archetypal Seyfert 2 galaxies, NGC 1068 and Mrk 3 
have both been the targets of recent X-ray grating obser- 
vations. The soft X-ray emission in both of these sources 
is extremely similar to the soft X-ray emission from NGC 
4151; it is extended on scales of hundreds of parsecs, shows 
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Figure 3. Best-fit of the ionization cone model to the RGS spectra of NGC 4151. The unbinned RGS 1 and RGS 2 spectra are shown 
in red and blue along with the best-fit model in green. The best-fit model parameters are given in Table 4. Unambiguously identified 
hydrogen-like and helium-like emission lines and RRC are labeled in a similar fashion to that in Figure 1. 



a strong correlation with the optical ionization cones iden- 
tified in HST O[III](A5007) images (Young et al. 2001; Sako 
et al. 2000) and the soft X-ray spectra are all dominated by 
similar emission lines and RRC. 

4.1 The comparison with the XMM-Newton RGS 
spectrum of NGC 1068 

The soft X-ray XMM-Newton RGS spectrum of NGC 4151 
is remarkably similar to that of NGC 1068. Both sources 
show strong hydrogen-like and helium-like line emission and 
RRC from carbon, nitrogen, oxygen, neon, magnesium and 
silicon. The resemblance between the spectra is apparent in 
the details of the spectra as well as on a more general level. 
For example, both sources show similar patterns of reso- 
nance, intercombination and forbidden line strengths for the 
helium-like triplets of nitrogen and oxygen and both show 
significant (but unresolved) iron L-shell emission. The prop- 
erties of the detected emission lines are remarkably alike in 
the two sources, both in terms of the lines present and the 
scale of their energy shifts. In both sources the prominent 
emission lines are blue-shifted by <1000 km s~^ with aver- 
age line shifts of 305 km s'^ and 375 km s~^ for NGC 4151 
and NGC 1068, respectively (Kinkhabwala et al. 2002). The 
line widths measured for NGC 1068 do not show the bi- 
modal distribution observed in NGC 4151. However, since 



only 7 lines have constrained line widths in the NGC 4151 
spectrum, compared to 23 lines with constrained widths in 
the NGC 1068 spectrum, caution should be strongly empha- 
sized in interpreting the validity of the distribution of line 
widths in NGC 4151. Interestingly the average line width 
found in the NGC 1068 data (~1100 km s"^ FWHM) is 
comparable to the average line width of the group of broader 
lines in the NGC 4151 spectrum (~1000 km s"^ FWHM). 
It is tempting, based upon a comparison of the helium-like 
triplet plasma diagnostic ratios (R and G; Section 3.3) to in- 
fer that the ionization cones in NGC 4151 have a somewhat 
lower electron temperature and a somewhat higher electron 
density than those in NGC 1068 (Kinkhabwala et al. 2002). 
Unfortunatly the situation may not be a clear as cut as this 
simple comparison suggests, primarily due to the unknown 
effect of photoexcitation on the diagnostic ratios (discussed 
previously in Seciton 3.3). The best-fit ionic column densi- 
ties for NGC 1068 were found to be similar to those observed 
m absorption for Seyfert 1 galaxies, implying either radially 
stratified ionization zones or, more likely, a broad distribu- 
tion of densities at any given radius. An inspection of the 
radial column densities fit to the RGS spectrum reveals that 
NGC 4151 appears to have less ionized material in its ion- 
ization cones, by a factor of '-^3.5 on average. The best-fit 
ionic electron temperatures, based on the widths of the RRC 
features, are also extremely similar between NGC 4151 and 
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NGC 1068 (although the comparatively large errors on these 
temperatures prevents any definite conclusions being drawn 
regarding this similarity) . The ionic column densities are de- 
pendent on a combination of factors including, cone length, 
density, ionization parameter and the filling factor of the 
material. Lower ionic column densities could be the result 
of lower ionizing flux or less material in the ionization cones, 
however a somewhat more plausible explanation, in line with 
the smaller R-values from the helium-like line triplets, could 
1)(> that the material in the ionization cones of NGC 4151 is 
rather more dense than the material in the ionization cones 
of NGC 1068. Since the recombination rate is dependent 
on the density of the material, this can result in a smaller 
ionized fraction for a similar ionizing continuum. 



4.2 The comparison with the Chandra HETG 
observation of Mrk 3 

The Chandra HETG observation of Mrk 3 was the first pub- 
lished high resolution soft X-ray spectrum of a Seyfert 2 
galaxy and showed unambiguously that the soft X-ray emis- 
sion in Mrk 3 is composed almost entirely of emission lines 
and RRC, with little or no continuum emission. The soft X- 
ray spectrum of Mrk 3 is, again, remarkably analogous to the 
soft X-ray spectrum of NGC 4151, though with somewhat 
less statistical significance due to the rapidly decreasing ef- 
fective area of the HETG instrument at wavelengths above 
~9A. In both sources there is strong evidence to suggest that 
the observed emission lines are produced in a photoionized 
and photoexcited plasma (Sako et al. 2000). The importance 
of photoexcitation in the spectrum of Mrk 3 is highlighted by 
the observed weak iron L-shell emission, rather than by the 
enhanced higher order resonance transitions, rejecting expla- 
nations from a recombination dominated or a collisionally 
dominated plasma. Despite this impressive spectrum, the 
Chandra observation of Mrk 3 only partially resolves two of 
the helium-like line triplets detected in the RGS spectrum of 
NGC 4151 (oxygen and neon), due to the limited number of 
photons in each feature, and as a result many of the proper- 
ties of the plasma are not well constrained. In particular, the 
resonance and intercombination lines of these triplets can- 
not be separated and as a result Sako et al. (2000) quote a 
different diagnostic line ratio {i.e. (r-|-i)/f) to the line ratios 
used here. 

At wavelengths shorter than ~9A the Chandra spec- 
trum shows evidence for a reflection dominated continuum 
along with clear detections of emission lines that are either 
extremely difficult or impossible to detect with the RGS 
(given its limited spectral bandpass), namely Fe XXVI, Fe 
XXV, Fe Kq, S XVI, S XV, S Kq, Si XIV, Si XIII, Si 
Ka, Mg XII and Mg XI (Sako et al. 2000). This plethora 
of lines highlights the complimentary capabilities of XMM- 
Newton and Chandra and strongly supports future obser- 
vations of bright, heavily obscured AGN with both Chan- 
dra (using the HETG) and XMM-Newton, to obtain the 
highest possible resolution X-ray spectra over the widest 
possible energy range. The line ratios calculated from the 
Chandra detections of the helium-like silicon and magne- 
sium line triplets (i.e. Si XIII and Mg XI) are all incon- 
sistent with the expected values for a collisionally ionized 
plasma and are marginally consistent with the predictions 
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Figure 4. Upper Panel: The 0.3-12 keV EPIC PN spectra 
of NGC 4151 modelled by a combination of the hard X-ray 
model developed in Schurch et al. (2003) and the 'ionization 
cone' soft X-ray model presented in Section 3.4. Lower Panel: 
The ratio of the data to the model prediction. 

from a photoionized plasma, similar to the conclusions from 
the triplet ratios for both NGC 4151 and NGC 1068. 



5 A BRIEF COMPARISON WITH THE 
XMM-Newton EPIC SPECTRA 

One of the strengths of XMM-Newton is that all the sci- 
ence instruments operate simultaneously during each ob- 
servation, allowing a direct comparison between the high 
resolution spectra from the RGS instruments and the lower 
resolution spectra from the EPIC CCD detectors. This com- 
parison is particularly useful in exposing features in each 
source spectrum that the other detectors struggle to detect. 
For example, broad features present in the 0.3-2.5 keV range 
arc not easily identified in RGS spectra, but will be more 
clearly observed by the EPIC instruments. Conversely, fea- 
tures that appear to be broad in the EPIC spectrum may 
be resolved into a set of individual narrow lines by the RGS 
instruments. Conclusions from the process of comparing the 
EPIC and RGS spectra are complicated by the need for 
good calibration and cross-calibration for each instrument 
involved. 

Figure 4 shows the full window mode EPIC PN spectra 
fit with a model that combines the best-fit of the hard X- 
ray model developed in Schurch et al. (2003) and the best-fit 
of the 'ionization cone' soft X-ray model presented in Sec- 
tion 3.4. In addition to the hard X-ray model and the ion- 
ization cone model, a silicon Kq line was included in the 
fitting. The energy of the line was fixed at 1.73 keV and the 
line width was fixed to be the same width as the iron Ka line 
width. The line flux was fixed to be 5% of the flux contained 
in the iron Ka line, in keeping with the ratio of the iron Ka 
to silicon Ka lino fluxes measured in the high resolution 
Chandra spectrum of Mrk 3. Initially the column densities 
of the hydrogen and helium-like atomic species were flxed to 
be the values measured by the RGS data. Key parameters 
in the hard X-ray model (e.g. ionization parameter, iron Ka 
lino flux etc; sec Schurch et al. 2003 for details) were free 
parameters in the fit, but no contributions from iron L-shell 
transitions were included in the fit. This highly restricted 
model does not provide a good fit to the EPIC spectrum. 
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with considerable residuals remaining below 2 koV. Relax- 
ing the constraints of the model by incorporating a contri- 
bution from iron L-shell emission, improves the fit consid- 
erably and results in similar values for the column densi- 
ties of iron ionization states from iron XVI-XXIII to those 
found with the RGS. Despite the improvement, the overall 
fit remains poor with considerable residuals still remaining 
across the soft X-ray region. Relaxing the constraints of the 
model more considerably by allowing the column densities 
of the hydrogen and heliurn-likc atomic species to be free 
parameters in the model fit, provides a reasonable fit to the 
broadband X-ray spectrum of NGC 4151 (x^=2709 for 2194 
d.o.f) however residuals remain in the fit below 2 kcV at a 
level of ~10% (Figure 4). The best-fit column densities for 
the hydrogen and liclium-like atomic species are typically 
within a factor of ~3 of the best-fit column densities de- 
termined from the RGS modelling. We do note quote these 
values here because the relatively poor spectral resolution 
of the EPIC instruments result in considerable degeneracy 
between the column densities measured for specific atomic 
species maJiing any specific values at best unreliable and at 
worst misleading. Despite this caveat, we note that the best- 
fit column densities of Nitrogen VIL Noon X and Silicon XIV 
measured in the EPIC data are considerably lower in com- 
parison with those measured in the RGS Data (factors of 3.5, 
40, and 16 respectively). The modelling is not significantly 
improved by the addition of either a soft X-ray power-law 
(representing a scattered fraction of the underlying contin- 
uum) or blackbody emission (to represent the canonical 'soft 
X-ray excess' often observed in Seyfert 1 galaxies). Consider- 
ably more detailed modelling of the RGS spectrum utilizing 
the additional data from the most recent XMM-Newton ob- 
servations of NGC 4151, coupled with improvements in the 
relative instrument calibrations, is required before any firm 
conclusions on the nature of the remaining fit residuals can 
be determined. 



6 SYNOPSIS 

We present a detailed analysis of the RGS spectrum of NGC 
4151 from a total of ~125 ks of XMM-Newton observations. 
The EPIC spectra from these observations are presented in 
Schurch et al. (2003). The analysis uses the high spectral 
resolution of the RGS instruments to characterize the nature 
and origin of the soft X-ray spectrum of NGC 4151. 

• The soft X-ray spectrum of NGC 4151 (Figure 1) is 
completely dominated by line emission and radiative recom- 
bination continua. Line emission, including higher order (/?- 
e) resonance transitions, is clearly detected from hydrogen- 
like and helium-like ionization states of neon, oxygen, ni- 
trogen and carbon. Much of the unidentified line emission 
below ~16A can be associated with a combination of iron L- 
shell omission and lines from hydrogen-like and helium-like 
silicon and magnesium. 

• The RGS spectra show no significant continuum emis- 
sion and no fluorescent line emission associated with neutral 
material. 

• The measured line centroid energies are consistent with 
a blueshift of between ^100-1000 km s~^, with respect to 
the systemic velocity of NGC 4151. The seven observed line 
widths that are constrained by the data tentatively suggest 



a birnodal distribution, with average line widths of '^lOOO 
km s"^ and ~200 km s"^ (FWHM). The measured line 
blueshifts are generally consistent with the outflow veloc- 
ities of the absorption lines observed in the UV spectrum 
of NGC 4151 (Kriss et al. 1995), implying that the UV and 
X-ray line features originate in the same material. 

• Plasma diagnostics, based on the ratios of the res- 
onance, intercombination and forbidden lines of the ob- 
served helium-like triplets, imply electron temperatures of 
~l-5x 10'* K and electron densities of ^^10* cm^'' in nitrogen 
and oxygen and ^10^" cm~^ from neon (Porquet & Dubau 
2000) We note again that the effects of photoexcitation may 
have a considerable effect on the values of electron tempera- 
ture and density determined by these diagnostic ratios (see 
Section 3.3). The wide range of the calculated electron den- 
sities is suggestive of either a 'clumpy' extended NLR or a 
region in which the material has grossly non-solar elemental 
abundances. 

• The soft X-ray spectrum of NGC 4151 is remarkably 
similar to the soft X-ray spectrum of NGC 1068 (Kinkhab- 
wala et al. 2002). The correspondence between the two 
source spectra is evident on the general scale from the 
strong emission lines and RRC associated with the same 
ionic species, which dominate the spectra of both sources. 
The resemblance of the NGC 4151 and NGC 1068 soft X-ray 
spectra is also apparent in the details of the spectra {i.e. the 
similar helium-like triplet ratios and RRC temperatures). 

• The 'ionization cone' model, developed by Kinkhabwala 
et al. (2002) to model the RGS spectra of NGC 1068, re- 
produces all of the hydrogen-like and helium-like emission 
features observed in the soft X-ray spectrum of NGC 4151 
in detail. The best-fit radial and transverse velocities of the 
material in the ionization cones agree well with the blueshift 
and velocity dispersion measured from the individual line 
centroid energies, and are similar to the radial and trans- 
verse velocities measured from the soft X-ray spectrum of 
NGC 1068. The individual ionic column densities arc some- 
what smaller than those measured for NGC 1068, which, 
in combination with the plasma diagnostic ratios from the 
helium-like triplets, suggests that the material in the ioniza- 
tion cones of NGC 4151 may be somewhat denser than the 
material in the ionization cones of NGC 1068. 
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